White dwarf cooling and large extra dimensions 
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Abstract 

Theories of fundamental interactions with large extra dimensions 
have recently become very popular. Astrophysical bounds from the Sun, 
red-giants and SN1987a have already been derived by other authors for 
the theory proposed by Arkani-Hamed, Dimopoulos and Dvali. In this 
paper we consider G117-B15A pulsating white dwarf (ZZ Ceti star) for 
which the secular rate at which the period of its fundamental mode 
increases has been accurately measured and claimed that this mode of 
G117-B15A is perhaps the most stable oscillation ever recorded in the 
optical band. Because an additional channel of energy loss (Kaluza- 
Klein gravitons) would speed up the cooling rate, one is able to use 
the aforementioned stability to derive a bound on theories with large 
extra dimensions. Within the framework of the theory with large extra 
dimensions proposed by Arkani-Hamed, Dimopoulos and Dvali we find 
the lower bound on string comapctification scale M s > 14.3 TeV/c 2 
which is more stringent than solar or red- giant bounds. 
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1 Introduction 



The interest in physical theories with extra spatial dimensions has recently ex- 
perienced considerable revival. Multidimensional theories with compact spa- 
tial dimensions having inverse radius at the order of the GUT scale were in- 
vestigated intensively in the 80-ties [p]]. Later developments in string theory 
initially due to Anoniadis |2| later supported by || suggested that it could 
be possible to have a much lower string or compactification scale. In par- 
ticular, it is has been conjectured 0] that compactification scale could be at 
the order of a TeV, corresponding to a weak-scale string theory. Such a low 
string scale is attractive from experimental perspective since the string state 
spectrum would now become accessible. In this class of theories, gravity is 
essentially n + 4 dimensional whereas all other physical fields are confined to 
4-dimensional brane. The relation between the Planck mass in 4 dimensions 
(Mpi = 1.2 10 19 GeV/c 2 ) and the string mass scale in 4 + n dimensions M s 
and the radius R of extradimensional space reads: 

R ~ \c) MTHT n (1) 

where Q n is the volume of the unit n-sphere. Present laboratory limits || give 
ITeV as lower bound on M s . Assuming that M s is of order of TeV (sustaining 
the hopes of experimental verification of multidimensionality of the world) one 
can immediately rule out n — 1 because in that case one would expect modified 
gravity at distances R ~ 10 15 cm which is not observed. 

Last century witnessed a successful application of standard physics in elu- 
cidating the properties of celestial bodies which consequently can be used as a 
source of - sometimes strong - bounds. The idea that astrophysical considera- 
tions can constrain "exotic" physics is not a new one. For example it has been 
implemented to constrain the axion mass ||. The main idea here is that if 
"exotic" physics to be tested predicts the existence of weakly interacting parti- 
cles which can be produced in stellar interiors such weakly interacting particles 
would serve as an additional source of energy loss in many ways influencing 
the course of stellar evolution. In the context of multidimensional physics 
excitation of Kaluza-Klein gravitons play the role of stellar extra coolant. 

There are three main sources of astrophysical bounds invoked in this con- 
text. First is the Sun which is a hydrogen burning main sequence star with 
radiative interior. Back reaction in response to increased energy loss results 
in rising the internal temperature and shrinking of the radiative interior (for 
details see ||). Because enhanced temperature increases the rates of nuclear 
reactions this would exacerbate solar neutrino problem. Moreover helioseismol- 
ogy provides an accurate estimate of the internal profile of isothermal sound 
speed squared M. The most accurate value of this quantity, which is pro- 
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portional to temperature, refers to the radius r = 0.2 R Q . This provides an 
effective mean to constrain exotic sources of energy loss. 

The second source of bounds comes from red giants - stellar evolutionary 
phase after the hydrogen has been exhausted and which lasts until the helium 
ignition takes place during the so called helium flash. Additional cooling would 
results in one of the following effects. First of all helium flash would not 
occur if cooling was too effective. If additional cooling is less effective (so 
that helium flash eventually takes place) the red giant tip on the Hertzsprung- 
Russel diagram would be located higher above the Horizontal Branch and the 
time spent on a Horizontal Branch would be shortened. All these effects could 
be tested on HR diagrams for globular clusters ||. 

The last and the most effective of traditional sources of bounds comes from 
SN1987a - more specifically from the observed duration of neutrino pulse. The 
pulse would be shortened had the nascent neutron star cooled down more 
rapidly than standard theory predicted. These considerations have been ap- 
plied, in the context of large extra-dimensions in works of Barger et al. [ IT ], 
Cullen and Perelstein [0 and Cassisi et al Jljj (see also fll3|). 



In this paper we consider another class of astrophysical objects for which 
the cooling rate is known from observations, namely the pulsating white dwarfs. 
We will focus our attention on ZZ Ceti star G117-B15A which was considered 
previously as a tool of testing fundamental physics (in the context of axion 
emission) || |(| . 



2 White dwarf cooling in the presence of large 
extra dimensions 

White dwarfs such like G117-B15A are final stages of evolution of stars with 
masses smaller than 8 M & . They have inner degenerate core composed of car- 
bon and oxygen which is an internal energy reservoir and a thin non-degenerate 
envelope made up of He and H shells. Having exhausted their nuclear fuel these 
stars can only contract and cool down. Some of them e.g. ZZ Ceti stars are 
pulsating and their pulsation period increases as the star cools down || [TO] . 



The rate of period increase P/P is proportional to temperature decrease rate 
f/T. 

At an energy scale much lower than the string scale, one can construct an 



effective theory of KK gravitons interacting with the standard model fields |15 



Barger et al. have calculated Kaluza-Klein graviton emissivities in five pro- 
cesses interesting from astrophysical perspective: photon-photon annihilation, 
electron - positron annihilation, Gravi-Compton-Primakoff scattering, Gravi- 
bremsstrahlung in a static electric field and nucleon-nucleon bremsstrahlung. 
Kaluza-Klein gravitons can couple to photons which leads to the first process, 
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similarly electron-positron pair can annihilate into Kaluza-Klein gravitons, in 
the third process scattering of photons by electrons may lead to Kaluza-Klein 
graviton emission via Compton or Primakoff processes (conversion of photons 
into Kaluza-Klein gravitons) - respective Feynaman diagrams can be found 
ITJ. First of the last two processes consists in bremsstrahlung emission of 
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Kaluza-Klein gravitons in static electric field of ions, the next one is similar - 
the bremsstrahlung emission of gravitons by nulceons (in electric field of nu- 
cleons). One can expect that first process becomes effective in hot stars where 
the photon density is large enough, the second one in stars with abundant 
electron-positron pairs (i.e. mainly the protoneutron or young neutron stars) 
etc. 

Because white dwarfs are dense and cool one can expect that dominant 
process of Kaluza-Klein graviton emission is gravi-bremsstrahlung of electrons. 
The specific (mass) emissivity estimated by Barger et al. |TT| is 



e = 5.86 10- 75 ^Vn 7 Z, 2 forn = 2 (2) 

pM 4 J 3 J V ) 

e = 9.74i0- 91 ^y>,Z 7 2 forn = 3 (3) 
pMf j 

where: T is the temperature of isothermal core, p is the density, n e and rij are 
the number densities of electrons and ions respectively. Total Kaluza-Klein 
graviton luminosity can be obtained as 

rM wn 

Lkk = edm (4) 

Jo 

Extensive theoretical studies of non-radial oscillations of white dwarfs had been 
carried out long time ago || . Appropriate theory of stellar oscillations consists 
in linearising the Poisson equation as well as equations of momentum, energy 
and mass conservation with respect to small non-radial perturbations. These 
perturbations (which in ZZ Ceti stars are the so called g-modes) can have os- 
cillatory behavior oc exp(— iat) where a 2 oc —A. The quantity A (which also 
illustrates the connection between oscillations and convective instability) is de- 
termined by thermodynamical properties of stellar matter: A = — 
where: p denotes density, r - radial coordinate, p is the pressure and I\ is 
the adiabatic index |TJ| . Now, for a zero-temperature degenerate electron gas 



A = meaning that no g-modes are supported. However, if non-zero thermal 
effects are taken into account one can show that A oc T 2 and consequently 
4 oc T i.e. the periods scale like 1/T where T is the core temperature. Con- 
sequently, the increase of the pulsation period can be calculated from the 
following formula: 

P T _ L 

P^'f- ~c v M WD T (5) 
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where L is the total luminosity, c v - average heat capacity and Mwd is the 
mass of the white dwarf. In the second part of the above formula (equality) 
the famous Mestel cooling law has been applied |L7|]. As mentioned above we 
assume that anomalous cooling rate observed in the form of pulsation period 
increase is due to Kaluza-Klein graviton emission. In such case we can write 

I 

LrK Fobs i /m 

- 1 (6) 
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where L 7 represents standard photon cooling, P b s represents the observed 
speed of pulsation period increase, P is an analogous quantity but without 
Kaluza-Klein cooling (theoretical). 

Since its discovery in 1976 |H| G117-B15A has been extensively studied. 
Regarding its variability the observed periods are 215.2, 271 and 304.4 s to- 



gether with higher harmonics and linear combinations thereof ||19|| . First esti- 
mates of the rate of increase of main pulsation period P = 215.2 s motivated 
Isern et al. M to discuss the effect of axion emission from this star. Very 
recently, with a much longer time interval of acquired data, Kepler et al. [^l| 
recalculated the rate of period increase and found significantly lower value of 
P = (2.3 ±1.4) x 1(T 15 s s -1 . Hence it has been claimed that the 215.2 s 
mode of G117-B15A is perhaps the most stable oscillation ever recorded in the 
optical band (with a stability compared to milisecond pulsars |p2[ ). 

This circumstance makes it possible to derive a bound on string mass scale 
M s . Previous upper limits on P would only allow to test whether such "pecu- 
liar" behavior (if true) could be a manifestation of large extra dimensions and 
whether it was consistent with other astrophysical bounds. 

The white dwarf pulsator G117-B15A has a mass of 0.59 M , effective 
temperature T e jf = 11620 K |33| and luminosity log(L/L Q ) = —2.8 [^4] (i.e. 



L 7 = 6.18 x 10 30 erg s^ 1 ). Typical model for such CO star predicts the central 
temperature T = 1.2 x 10 7 K and matter density p = 0.97 x 10 6 g cm -3 . 
We assume that mean molecular weight per electron is fi e ~ 2. This allows 
us to estimate electron number density n e . Following detailed calculations 
by Salaris |25[ chemical composition of the white dwarf core has been taken 



as 83% by mass of oxygen and remaining 17% of carbon. Then by virtue of 
electric neutrality of the star one is able to estimate rij. Now we can estimate 
the power radiated away from the white dwarf core in the form of Kaluza-Klein 
gravitons Lkk- In order to derive the lower bound on M s one can take (in 
the role of P a b s in (|)) an upper 2a limit of P equal to 5.1 x 10~ 15 s s^ 1 |2~T| . 



Following [E3 we assume Pq = 3.9 x 10 15 s s This means that recently 



established secular stability of the fundamental oscillation mode of G117-B15A 
implies Lkk < 0.308 L 7 which translates to: 

M s > 14.3 TeV/c 2 forn = 2 
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Respective graviton emission rates for n = 3 (and greater) theories turn 
out to be negligible, hence we do not quote the resulting numbers. 

Our bound has been derived within a very simple law relating the observed 
rate of the pulsational period and the the rate of the change of the period 
given by models when Kaluza-Klein gravitons are considered - essentially that 
proposed by Isern et al. ||. The most self-consistent approach would be to 
run a set of evolutionary models of white-dwarfs and compare their predictions 
concerning P with the observed value P bs- In the case of axion emission this 
has been done by Corsico et al. in and the result obtained that way turned 
out to be in a very good agreement with the one obtained by applying the above 
mentioned much simpler method. The reason for this could be understood by 
very simple argument that the axion (or Kaluza-Klein gravitons as in our 
case) emissivity is dominated by the bremsstrahlung process in isothermal 
degenerate core. Based on the evolutionary white dwarf models the authors 
of [^2| pointed out to a very important circumstance that additional cooling 
from the "exotic" physics affects the temperature profile in the innermost 
degenerate parts of the star. The structure of the outer (partially degenerate) 
layers depends on the temperature profile. However this profile is to great 
extent determined by the effective temperature (fixed by observations). In 
consequence the periods of oscillations (which depend on the structure of the 
star) remain almost unchanged when additional cooling (which in any case is a 
small correction to the energy budget) is present. On the other hand the value 
of P is sensitive to additional cooling. This justifies the assumption that one 
can first identify the structure of the fiducial model (taken by us as that cited 
in P2J) without Kaluza-Klein graviton emission and then incorporate graviton 
emission in discussing the rate of secular changes of the period. 



3 Conclusions 

It is interesting to compare our result with existing astrophysical constraints 
on the string compactification scale within the framework of the theory pro- 
posed by Arkani-Hamed, Dimopoulos and Dvali H. Helioseismological and 



red-giant type considerations were performed by Cassisi et al. [13 1. They 
calculated detailed solar models taking into account energy loss in Kaluza- 
Klein gravitons explicitly in their code and obtained the lower bound for M s 
equal to 0.3 TeV/c 2 . More stringent limit was derived from simulating the 
globular cluster Hertzsprung-Russel diagram [|14]] (implementing Kaluza-Klein 
graviton emissivity into FRANEC evolutionary code). By virtue of comparing 
predicted luminosity of RGB tip with observations Cassisi et al. obtained a 
"red-giant" bound for M s to be 3 — 4 TeV/c 2 . Our estimate of the string 
energy scale implied by recently reported stability of ZZ Ceti star G117-B15A 
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equal to 14.3 TeV/c 2 is much stronger than above mentioned stellar evolu- 
tionary bounds. 

Among existing astrophysical bounds on Kaluza-Klein theories with large 
extra dimensions only the supernova constraints are more restrictive. They 
demand M s > 30 — 130 TeV/ c 2 \TL, |T^| and are based on a different mechanism 



of Kaluza-Klein graviton emission - the nucleon-nucleon bremsstrahlung. On 
the other hand the lesson learned in testing the physics of axions shown that 
first straightforward supernova bounds were reduced by an order of magnitude 
when more accurate nuclear physics was employed f27j. Quite recently, the 
paper with improved calculations of nucleon-nucleon bremsstrahlung appeared 



13 1 suggesting that former SN1987a bounds should be lowered by a factor of 



0.65. It should also be noted that cosmological considerations [f28j apparently 
provide much more stringent bounds than the supernova SN1987a. 

The bound derived in this paper is based on white dwarf cooling. The 
physics underlying this process is very simple hence one can expect that the 
result is robust (within the framework of the theory proposed by Arkani- 
Hamed, Dimopoulos and Dvali |4]]). Kaluza-Klein graviton emissivity at the 
relevant densities and temperatures of white dwarfs is dominated by the gravi- 
bremsstrahlung process taking place in the degenerate and isothermal core. 
The mass of the core (essentially equal to the mass of the star, since the outer 
helium layer comprises less than 0.01 M^b, its temperature and chemical 
composition can be reliably estimated by fitting evolutionary models to obser- 
vational characteristics such like effective temperature or oscillation periods. 
The cooling rate of B117-G15A pulsating star is constrained by measurements 
which are performed with great accuracy. The most recent determination of 
the secular rate of change of the period |21| took into account all the peri- 
odicities and the error bars reported therein can be considered as safe. It 
can be argued [22] that remaining uncertainty - mostly from mode identifica- 
tion procedure and the precise physical characteristics (mass or temperature) 
is of order of 1 x 10~ n s s -1 . Other effects like the contribution of the 
proper motion p(| |21| or the rate of reaction 12 C(a,7) 16 which determines 
the stratification of the core at final stage of the asymptotic giant branch, both 
contribute an order of magnitude smaller value to the final uncertainty. 

The detailed discussion of using a pulsational code coupled to evolutionary 
code aimed at constraining an axion mass by observed stability of the funda- 
mental mode of G117-B15A can be found in a recent paper by Corsico et al. 
22 1 . One of the conclusions formulated in p2| was that the bounds on axion 



mass derived from simple estimates like performed in the present paper (in a 
different context) or in || are in good agreement with evolutionary calcula- 
tions. Although axion emissivity has different temperature dependence than 
that of Kaluza-Klein gravitons, one can expect the same for gravitons. Hence 
the B117-G15A pulsator remains an important tool for testing fundamental 
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physics. 
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